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Abstract 

A search for the Standard Model Higgs boson in the two photon decay channel is reported, using 1.08 fbr 1 of 
proton-proton collision data at a centre-of-mass energy of 7 TeV recorded by the ATLAS detector. No significant 
excess is observed in the investigated mass range of 1 10-150 GeV. Upper limits on the cross-section times branching 
ratio of between 2.0 and 5.8 times the Standard Model prediction are derived for this mass range. 



1. Introduction 

The search for the Standard Model Higgs boson [Q]- 
[3] is one of the key goals of the Large Hadron Collider 
(LHC) at CERN. The allowed Higgs boson mass (/«#) is 
constrained at the 95% confidence level by a lower limit 
of 1 14.4 GeV from the LEP experiments [4] and an ex- 
cluded region between 156 and 177 GeV from the Teva- 
tron experiments |5j |6J] . First results have been reported 
by the ATLAS experiment in a variety of channels 10] 
and the CMS experiment fit] using the data recorded 
in 2010, which correspond to an integrated luminosity 
about thirty times smaller than the 201 1 dataset used in 
this analysis. 

In the low mass range, from the LEP limit to mg 
» 140 GeV, one of the most promising search chan- 
nels at the LHC is the rare decay of the Higgs boson 
into a pair of photons. Despite the low branching ra- 
tio (« 0.2%), this channel provides good experimental 
sensitivity in the mass region below 150 GeV. The re- 
sults presented in this paper are based on proton-proton 
collision data taken at -\fs = 7 TeV by the ATLAS ex- 
periment between April and June 201 1 . 

The data analysis proceeds by selecting photon pairs 
with tight identification and isolation cuts to minimize 
backgrounds other than direct diphoton production. A 
narrow peak in the reconstructed invariant mass distri- 
bution is searched for over a large, smooth background 
whose normalisation and shape are left free in a maxi- 
mum likelihood fit. To increase the sensitivity, the sam- 
ple is divided into five categories based on the presence 
of photon conversions and on the photon impact point 
on the calorimeter, with different invariant mass reso- 
lutions and signal-to-background ratios for the different 



categories. 

The results of the fit are compared to the prediction 
from the Standard Model using the Higgs boson pro- 
duction cross-section and branching ratio from Ref. |^]. 
Limits on the production cross-section relative to the 
Standard Model value are then derived as a function 
of the hypothesised Higgs boson mass. Although with 
the current dataset the analysis is not yet sensitive to 
the predicted rate for a Standard Model Higgs boson, 
the limits on the yield in this decay channel improve 
on those obtained in the same channel by the Tevatron 
experiments ifiol- 121. and are sensitive to possible en- 
hancements in the Higgs boson production and decay 
rate compared with the Standard Model expectations. 

2. Experimental setup and data set 

The ATLAS detector is described in detail in 
Ref. lfl3tl . The main subdetectors relevant to this analy- 
sis are the calorimeter, in particular its electromagnetic 
section, and the inner tracking detector. 

The electromagnetic calorimeter is a lead-liquid ar- 
gon sampling calorimeter with accordion geometry. It 
is divided into a barrel section covering the pseudora- 
piditjQ region \r/\ < 1.45 and two end-cap sections cov- 
ering the pseudorapidity region 1.375 < \q\ < 3.2. It has 



ATLAS uses a right-handed coordinate system with its origin at 
the nominal interaction point (IP) in the centre of the detector and the 
z-axis along the beam pipe. The A-axis points from the IP to the centre 
of the LHC ring, and the y axis points upward. Cylindrical coordinates 
(r, <f>) are used in the transverse plane, <p being the azimuthal angle 
around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 9 as r\ = - In tan(#/2). 
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three longitudinal layers. The first one, with a thickness 
between 3 and 5 radiation lengths, has a high granu- 
larity in rj (between 0.003 and 0.006 depending on r], 
with the exception of the regions 1.4 < \rj\ < 1.5 and 
\rj\ > 2.4), sufficient to provide discrimination between 
single photon showers and two photons from a jt° decay. 
The second layer has a thickness of around 17 radiation 
lengths and a granularity of 0.025 x 0.025 in 77 x (p. A 
third layer, with a thickness varying between 4 and 15 
radiation lengths, is used to correct for leakage beyond 
the calorimeter for high energy showers. In front of the 
calorimeter, a thin presampler layer, covering \rj\ < 1.8, 
is used to correct for fluctuations in upstream energy 
losses. The sampling term a of the energy resolution, 
<t(E)/E * a I V£ (GeV), varies between 9% and 14% as 



a function of \q\ for unconverted photons [ 14]. It reaches 
up to 20% for converted photons near I77I of 1.3 where 
the upstream material effect is the largest. The sam- 
pling term is the largest contribution to the resolution 
up to about 100 GeV, where the constant term starts to 
dominate. After 0.17 ftr 1 of data were accumulated, 
some calorimeter cells could not be read out. The af- 
fected region size is A/7 x A(px 1 .5 x0.2 in the barrel elec- 
tromagnetic calorimeter, resulting in an acceptance loss 
for diphoton candidates of about 3%. A hadronic sam- 
pling calorimeter is located behind the electromagnetic 
calorimeter. It is made of steel and scintillating tiles in 
the barrel section, and of copper and liquid argon in the 
end-cap. 

The inner detector consists of three subsystems: at 
small radial distance R from the beam axis (5 < R < 15 
cm), pixel silicon detectors are arranged in three cylin- 
drical layers in the barrel and in three disks in each end- 
cap; at intermediate radii (30 < R < 56 cm), double lay- 
ers of single-sided silicon microstrip detectors are used, 
organised in four cylindrical layers in the barrel and nine 
disks in each end-cap; at larger radii (56 < R < 107 cm), 
a straw tracker with transition radiation capabilities is 
used. These three systems are immersed in a 2 T axial 
magnetic field. The silicon pixel and strip subsystems 
cover the range \r/\ < 2.5, while the transition radiation 
tracker acceptance is limited to the range \rj\ < 2.0. The 
inner detector allows reconstruction of secondary ver- 
tices, in particular of photon conversions occurring in 
the inner detector material up to a radius of w 80 cm. 

The total amount of material in front of the first active 
layer of the electromagnetic calorimeter (including that 
in the presampler) varies between 2.5 and 6 radiation 
lengths as a function of pseudorapidity, excluding the 
transition region (1.37 < \t]\ < 1.52) between the barrel 
and the end-caps. 

Data used in this analysis were selected using a di- 



photon trigger with a 20 GeV transverse energy thresh- 
old on each photon. At the first trigger level, which uses 
reduced granularity, two clusters with transverse ener- 
gies above 14 GeV are required in the electromagnetic 
calorimeter. At the higher trigger levels, loose photon 
identification cuts are applied using the full calorimeter 
granularity. This trigger has an efficiency greater than 
99% for the signal after the final event selection. 

In these data, the instantaneous luminosity varies be- 



tween sil0 32 cm _/! s _1 and =sl0 :i1 cm _2 s _l with a bunch 
spacing of 50 ns. The average number of collisions 
per bunch crossing is around 6. Collisions in the same 
bunch crossing as the signal (in-time pileup) or in other 
bunch crossings within the detector sensitive time (out- 
of-time pileup) influence the event reconstruction. The 
inner detector is only sensitive to in-time pileup while 
the electromagnetic calorimeter is sensitive to pileup 
within a as 450 ns time window. 

The application of beam, detector, and data-quality 
requirements to the recorded data results in a data sam- 
ple corresponding to a total integrated luminosity of 
(1.08+0.04) ft- 1 |l5]. 



3. Simulated samples 

The Higgs boson signal from the dominant gluon 
fusion production process (corresponding to 86% 
of the production cross-section for a Higgs bo- 
son with a mass of 120 GeV) is generated with 
POWHEG id. MC@NLO (lj is used as a cross- 
check. POWHEG d is also used to generate the signal 
events from the sub-leading vector boson fusion process 
(7% of the cross-section at 120 GeV). For the other pro- 
duction modes, namely associated production with a W 
or Z boson or a it pair, PYTHIA IU9I1 is used. 

The predicted signal is normalised using NNLO 
cross-sections for the gluon fusion process ll20H24ll . the 
vector boson fusion process 02511 . the associated produc- 
tion with a W or Z boson 12611 and NLO cross-section 
for the associated production with a tt pair 0. The 
NLO electroweak corrections are applied to the gluon 
fusion [28, 29|, vector boson fusion |30[ 31], and the 



associated production with a W or Z boson ll32Tl pro- 
cesses. The uncertainty on the theoretical cross-section 
is estimated Jgt] to be ^5%, mostly due to the renormal- 
isation and factorisation scale variations and the uncer- 
tainties in the parton distribution functions [33-36]. The 
Higgs boson decay branching fractions are taken from 
Refs. HIH. The uncertainty on the branching ratio 
to two photons is negligible compared with the cross- 
section uncertainty. 
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Signal events are generated in steps of 5 GeV for 
Higgs boson masses in the range of 110-150 GeV. 
PYTHIA and ALPGEN Q have been chosen to gen- 
erate the background samples, which are, however, only 
used for cross-checks and not to extract the final results. 

All Monte Carlo (MC) samples are processed through 
a complete simulation of the ATLAS detector 1I39I1 us- 
ing the GEANT4 programme [40]. Pileup effects are 
simulated by overlaying each MC event with a variable 
number of MC inelastic pp collisions, taking into ac- 
count both in-time and out-of-time pileup and the LHC 
bunch train structure. MC events are weighted to have 
the same distribution of average number of interactions 
per bunch crossing as in the data. 

4. Photon reconstruction, event selection and back- 
grounds 

4.1. Photon reconstruction 

Photon reconstruction is seeded by energy clusters 
in the electromagnetic calorimeter with transverse en- 
ergies exceeding 2.5 GeV in projective towers of size 
0.075x0. 125 in 77 x made from the presampler and the 
three electromagnetic calorimeter layers. These energy 
clusters are then matched to tracks that are reconstructed 
in the inner detector and extrapolated to the calorimeter. 
Clusters without matching tracks are classified as un- 
converted photon candidates. Clusters matched to either 
pairs of tracks which are consistent with the hypothesis 
of a photon conversion or single tracks without hits in 
the pixel layer nearest to the beam pipe are considered 
as converted photon candidates. The photon reconstruc- 
tion efficiency is «98%. 

The energy measurement is made in the electromag- 
netic calorimeter using a cluster size which depends 
on the photon classification. In the barrel, a size of 
0.075x0.125 in 77 x is used for unconverted photons 
and 0.075x0.175 for converted photon candidates, to 
account for the larger spread of the shower in <p for con- 
verted photons due to the magnetic field. In the end- 
cap, a cluster size of 0.125x0.125 is used for all can- 
didates. A dedicated energy calibration [ 14] is applied 
to account for upstream energy losses, lateral leakage 
and longitudinal leakage, separately for converted and 
unconverted photon candidates. 

The final energy calibration is determined from Z — > 
ee decays, resulting in 77-dependent correction factors of 
the order of ±1%. After this calibration procedure, the 
constant term in the energy resolution is estimated to be 
l.l*Qg% in the barrel region and 1.8^ ?% in the end-cap 



region [41]. The energy resolution in the simulation is 
adjusted to match these values. 



Photon identification is based on the lateral and lon- 
gitudinal energy profiles of the shower in the calorime- 
ter Il42ll . The photon candidate is required to deposit 
only a small fraction of its energy in the hadronic 
calorimeter. The transverse shower shape in the sec- 
ond layer of the electromagnetic calorimeter needs to 
be consistent with that expected for a single electromag- 
netic shower. Finally, the high granularity first layer is 
used to discriminate single photons from overlapping 
photon pairs from neutral meson decays produced in jet 
fragmentation, which are the main background source. 
Based on these criteria, a set of tight identification cuts, 
different for converted and unconverted candidates, is 
applied. 

To take into account small differences in shower 
shapes between data and simulation, the shape variables 
are shifted in the simulation before the identification 
cuts are applied. The photon identification efficiency 
ranges typically from 75% to 90% for transverse ener- 
gies between 25 and 100 GeV. 

To increase the background rejection, an isolation cut 
is applied. The isolation variable [42] is computed by 
summing the transverse energy in calorimeter cells in a 
cone of radius 0.4 in the 77 x <p space around the pho- 
ton candidate. Cells in the electromagnetic calorimeter 
within 0.125x0.175 from the shower barycentre are ex- 
cluded from the sum. The small photon energy leak- 
age outside the excluded cells is evaluated as a func- 
tion of the transverse energy in simulated samples and 
is subtracted from the isolation variable. To reduce the 
effect from the underlying event and pileup, the isola- 
tion is further corrected using a method suggested in 
Ref. 114 311 : for each of the two different pseudorapidity 
regions \t]\ < 1.5 and 1.5 < I77I < 3.0, low energy jets 
are used to compute an "ambient" energy density, which 
is then multiplied by the area of the isolation cone and 
subtracted from the isolation energy. 

In the following, photon candidates having isolation 
transverse energies lower than 5 GeV are considered as 
isolated. The isolation cut efficiency is checked in data 
using a control sample of Z — > ee events. The per-event 
efficiency of requiring both electrons to be isolated is 
found to be 3% lower in the data than in the simulated 
samples. In the MC, the isolation cut efficiency is found 
to be the same for Z — > ee and H — > yy events (a;93%). 
The number of events predicted by the simulation after 
the isolation cut is therefore reduced by 3%. 

4.2. Event selection 

Two photon candidates are required to pass tight 
identification criteria, to be isolated, and to be within the 
region I77I <2.37, excluding 1.37< I77I <1.52, where the 
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first calorimeter layer has high granularity. The high- 
est and second highest photon transverse energies are 
required be above 40 and 25 GeV respectively. Both 
photons must be clear of problematic regions in the 
calorimeter. As the goal is to investigate Higgs boson 
mass hypotheses between 110 and 150 GeV, the invari- 
ant mass of the photon pair is required to be within 100- 
160 GeV. After these cuts 5063 events remain in the se- 
lected data sample. 

The acceptance of the kinematic cuts, as estimated 
with generated photons in the MC signal samples, is 
60% for the dominant gluon fusion process for a mass 
of 120 GeV. The overall event selection efficiency, tak- 
ing into account both kinematic cut acceptance and re- 
construction and identification efficiencies, is 39%. The 
event selection efficiency is slightly larger in the vector 
boson fusion process. It is somewhat smaller in the as- 
sociated production mode. It increases with the Higgs 
boson mass from 34% at 1 10 GeV to 43% at 150 GeV. 

To enhance the sensitivity of the analysis, the data 
sample is split in five categories, with different invari- 
ant mass resolutions and different signal-to-background 
ratios: 

• Unconverted central (8% of the candidates): Both 
photons are unconverted and in the central part of 
the barrel calorimeter Qt]\ <0.75). This is the cat- 
egory with the best invariant mass resolution and 
the best signal-to-background ratio; 

• Unconverted rest (28% of the candidates): Both 
photons are unconverted and at least one pho- 
ton does not lie in the central part of the barrel 
calorimeter; 

• Converted central (7% of the candidates): At least 
one photon is converted and both photons are in the 
central part of the barrel calorimeter; 

• Converted transition (16% of the candidates): At 
least one photon is converted and at least one pho- 
ton is near the transition between barrel and end- 
cap calorimeter (1.3< \rj\ < 1 .75). Given the larger 
amount of material in this region, the energy reso- 
lution, in particular for converted photons, can be 
significantly degraded; 

• Converted rest (41% of the candidates): All other 
events with at least one converted photon. 

Table Q] shows the cross-section times branching ra- 
tio to two photons, the expected total and per category 
signal yields for 1.08 fbr 1 for different Standard Model 
Higgs boson mass hypotheses. Using these categories 



improves the signal sensitivity of the analysis by around 
15% for a 120 GeV Higgs boson mass compared with a 
fully inclusive analysis. 



Table 1 : Cross-section times branching ratio and expected numbers of 
signal events after all cuts (total and per category), for various Higgs 
boson masses and for an integrated luminosity of 1.08 fb _1 . 



mji [GeV] 


110 


120 


130 


140 


150 


crxBR [fb] 


45 


43 


37 


27 


16 


Signal yield 


17.0 


17.6 


15.8 


12.1 


7.7 


Unconverted central 


2.6 


2.6 


2.3 


1.7 


1.1 


Unconverted rest 


4.6 


4.7 


4.2 


3.4 


2.1 


Converted central 


2.0 


2.0 


1.7 


1.3 


0.8 


Converted transition 


2.3 


2.2 


2.1 


1.5 


1.0 


Converted rest 


5.6 


6.0 


5.6 


4.2 


2.7 



4.3. Invariant mass reconstruction 

In addition to the energies, the angle between the pho- 
tons is needed for the computation of the diphoton in- 
variant mass. This angle is determined from the interac- 
tion vertex position and the photon impact points in the 
calorimeter. The resolution of the angle measurement is 
dominated by the reconstruction of the primary vertex 
z position. The RMS vertex spread in the z direction is 
«5.5 cm, and a more accurate event-by-event estimate 
is performed to reduce the impact on the invariant mass 
resolution. Given the non-negligible level of pileup in 
the 2011 data, the determination of the vertex position 
is based only on the photon candidates, without relying 
on other charged tracks in the event. For converted pho- 
tons with tracks having a precise measurement in the z 
direction, the vertex position is estimated from the in- 
tercept of the line joining the reconstructed conversion 
position and the calorimeter impact point with the beam 
line. For all other photons, the vertex position is es- 
timated from the shower position measurements in the 
first and second layers of the electromagnetic calorime- 
ter, which can be used to calculate the photon direction. 
Finally, the vertex positions from both photons are com- 
bined taking also into account the average beam spot 
position in z. When both photons are unconverted, the 
typical vertex position resolution is a: 1.6 cm in z. The 
resolution is better in events with converted photons. 
The resulting impact of the angle measurement on the 
invariant mass resolution is negligible compared to the 
contribution from the photon energy resolution. 

Figure Q] shows the invariant mass distribution for 
simulated Higgs boson events with mass 120 GeV. The 
mass resolution for the signal is modelled by the sum of 
a Crystal Ball function 114411 (for the bulk of the events 
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which have a narrow Gaussian spectrum in the peak re- 
gion and tails toward lower reconstructed mass) and a 
Gaussian distribution with a wide sigma (to model the 
far outliers in the distribution). The Crystal Ball func- 
tion is defined as: 

| e-' 2/2 if t > -a C B 

i (HciLfcB . e -« 2 CB l^ .(!!£§__ t) -n C B otherwise 

where t = (m yy - /J-cb)/o~cb, N is a normalisation 
parameter, pcb is the peak of the narrow Gaussian 
distribution, <x CB represents the Gaussian resolution for 
the core component, and uqb and oqb parametrise the 
non-Gaussian tail. 

The core component of the mass resolution, cr CB , 
ranges from 1 .4 GeV in the "Unconverted central" cat- 
egory to 2.1 GeV in the "Converted transition" cate- 
gory. The non-Gaussian contributions to the mass reso- 
lution arise mostly from converted photons with at least 
one electron losing a significant fraction of its energy 
through bremsstrahlung in the inner detector material. 




<% [GeV] 



Figure 1 : Distribution of the reconstructed diphoton invariant mass of 
a simulated 120 GeV mass Higgs boson signal, for all categories to- 
gether. The line shows the fit of the mass resolution using the function 
described in the text. The core component of the mass resolution is 
1.7 GeV 



purity of the sample are described in more detail. Al- 
though the final result does not rely on it, a quantitative 
understanding of the sample composition is an impor- 
tant cross-check of the diphoton selection procedure. 

A method based on the use of control regions for two 
discriminating variables is applied to measure the con- 
tributions of fake photon background directly from the 
data. This method exploits relaxed isolation and pho- 
ton identification cuts to estimate the fake components, 
by relying on the fact that the rejections from these two 
cuts are almost independent. It is a generalisation of the 
method used in Ref. [42]. The Drell-Yan background is 
estimated by measuring the probability for an electron 
to be reconstructed as a photon candidate with Z events 
and applying this probability to the observed yield of 
Drell-Yan events at high mass. 

The number of diphoton events in the 100-160 GeV 
mass range is found to be 3650+100+290, where the 
first uncertainty is statistical and the second is system- 
atic. The systematic uncertainty arises from the defini- 
tion of the relaxed identification control region, the pos- 
sible correlations between isolation and identification 
variables, and the fraction of real photons leaking into 
the background control regions. The extracted yields of 
photon-jet and dijet are 1 1 10+60+270 and 220+20+130 
events respectively. The Drell-Yan background, which 
is most prominent in the categories with at least one 
converted photon, is estimated to be 86+1 + 14 events 
in the mass range of 100-160 GeV. 

Figure [2] shows the extracted components of the 
diphoton, photon-jet, dijet and Drell-Yan processes. 
The purity of the sample (fraction of diphoton events) 
is about 72%. The measurement of the purity has also 
been made separately in each category, and ranges from 
69% to 83%. 

Other methods have been used to cross-check the pu- 
rity estimate, in particular using template fits of the pho- 
ton isolation distribution, where both signal and back- 
ground templates are derived from data. The results are 
in agreement with the results quoted here. 



4.4. Sample composition 

The main background components are the diphoton 
production, the photon-jet production with one fake 
photon from jets fragmenting into a high energy n°, 
the dijet production with two fake photons, and Drell- 
Yan events where both electrons are misidentified as 
photons. A measurement of the diphoton production 
cross-section with 2010 ATLAS data can be found in 
Ref. ll45ll . where the techniques used to estimate the 



5. Systematic uncertainties 

Experimental systematic uncertainties affecting the 
extraction of the signal from the diphoton invariant mass 
distribution related to the modelling of the signal can be 
classified in two types: uncertainties affecting the pre- 
dicted yield and uncertainties affecting the modelling of 
the mass resolution. 

The uncertainties on the event yield are the following: 
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Figure 2: Diphoton, photon-jet, dijet and Drell-Yan contributions to 
the diphoton candidate invariant mass distribution, as obtained from 
a data-driven method. The various components are stacked on top of 
each other. The error bars correspond to the uncertainties on each 
component separately. 



The uncertainty from the photon reconstruction 
and identification efficiency amounts to ±11% per 
event. It is estimated from data and MC differences 
in shower shape variables, the impact of additional 
material in front of the calorimeter and the impact 
of pileup on the photon shower shape variables. 

The uncertainty on the isolation cut efficiency is 
taken as the difference between data and MC found 
in Z — » ee decays and amounts to +3% per event. 

The uncertainty on the photon trigger efficiency is 
±1%. It comes from the uncertainty in the mea- 
surement of the trigger efficiency for diphoton can- 
didates using control triggers and from possible 
differences between the trigger efficiency for pho- 
tons from Higgs boson decays and all diphoton 
candidates. 

The uncertainty on the kinematic cut acceptance 
from the modelling of the Higgs boson trans- 
verse momentum distribution is investigated with 
HQT [46] and RESBOS 0H, which account 
for all-orders soft-gluon resummation up to NNLL 
accuracy. The resulting uncertainty is found to be 
at the level of + 1%. 



The luminosity uncertainty is 3.7% 11511 . 



The total uncertainty on the expected signal event yield 
is +12%. 

The uncertainties on the invariant mass resolution are 
the following: 



• The uncertainty on the cluster energy resolution 
comes from the uncertainty on the sampling term, 
estimated to be 10%, and from the uncertainty on 
the constant term, which is estimated using Z — > ee 
decays. Both uncertainties are taken into account 
with their proper correlation from the Z control 
sample constraint. The uncertainty on the cluster 
energy resolution amounts to a +12% relative un- 
certainty on the diphoton invariant mass resolution. 

• The uncertainty on the photon energy calibration 
arising from the extrapolation of the electron en- 
ergy scale calibration is estimated from MC stud- 
ies. The difference between the photon and the 
electron response in the calorimeter comes from 
the material in front of the active part of the 
calorimeter. The uncertainty is estimated using 
simulations with a different amount of material in 
front of the calorimeter and is found to be ±6% on 
the mass resolution. 

• The contribution of pileup fluctuations to the clus- 
ter energy measurement is checked using random 
clusters in randomly triggered bunch crossings, 
with a frequency corresponding to expectations 
from the instantaneous luminosity. The relative un- 
certainty on the mass resolution is found to be less 
than 3%. 

• The uncertainty on the resolution of the photon an- 
gle measurement is studied with Z — > ee decays. 
The calorimeter-based direction measurement is 
compared with the much more precise track-based 
direction measurement. In the barrel calorime- 
ter, the resolution measured in data agrees well 
with the one predicted by the simulation. In the 
end-cap region, the resolution measured in data is 
«20% worse than in the simulation. The impact of 
this difference is a 1% relative uncertainty on the 
diphoton mass resolution. 

The total relative uncertainty on the diphoton invariant 
mass resolution is thus +14%. This systematic uncer- 
tainty is applied to both the Crystal Ball and the wide 
Gaussian resolution parameters. 

These systematic uncertainties are taken as fully cor- 
related between the different categories. The impact 
of uncorrelated systematic uncertainties in the different 
categories and migration between categories has been 
investigated and found to be negligible. 

The background is modelled by an exponentially 
falling invariant mass distribution. Systematic uncer- 
tainties arise from possible deviations of the back- 
ground distribution from this assumed shape. This has 
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been estimated by checking how accurately the cho- 
sen model fits different predicted diphoton mass dis- 
tributions ll49L l50ll and comparing different functional 
forms for the background model. The resulting uncer- 
tainty is between +5 events at 110 GeV and ±3 events 
at 150 GeV for a Higgs boson mass signal region about 
4 GeV wide. 



6. Results 

The data are compared to background and signal- 
plus-background hypotheses using a profile likelihood 
test statistic as described in Refs. fx 51]. The expo- 
nentially falling invariant mass distribution used for the 
background model is determined by two nuisance pa- 
rameters per category (the normalisation and the expo- 
nential negative slope), which are left free in the un- 
binned fit. The signal is modelled with the mass resolu- 
tion functions described above, one per category, fixing 
the fraction of events in each category to the MC pre- 
dictions. Table [2] summarises the measured fractions of 
background events in each category, the predicted frac- 
tions of signal events and the predicted core gaussian 
signal mass resolutions for a Higgs boson mass hypoth- 
esis of 120 GeV. 



Table 2: Fractions of background (ft), predicted signal (/,) and core 
gaussian mass resolution (<x) in each category for a Higgs boson mass 
hypothesis of 120 GeV. The fractions are normalised to the total yield 
summed over categories. 



category 


A 


/; 


cr(GeV) 


Unconverted central 


1% 


15% 


1.4 


Unconverted rest 


29% 


27% 


1.6 


Converted central 


8% 


11% 


1.5 


Converted transition 


16% 


13% 


2.1 


Converted Rest 


40% 


34% 


1.8 



The fitted parameters for the signal are thus the over- 
all signal strength relative to the Standard Model pre- 
diction and the nuisance parameters on the predicted 
event yield and mass resolution which have Gaussian 
constraints in the fit. The uncertainty on the predicted 
event yield includes both the experimental systematic 
uncertainties described in section [5] and the uncertainty 
on the predicted cross-section described in section [3] 
The systematic uncertainty on the background shape is 
included as another nuisance parameter with a Gaussian 
constraint in the fit. From this fit, the best estimate of the 
signal yield is extracted, as well as the likelihood ratio 
(profile likelihood) between any assumed signal yield 
(leaving the nuisance parameters free to maximise the 



likelihood) and the best estimate. The fit is performed 
in 1 GeV steps for the Higgs boson mass hypothesis, 
which is significantly smaller than the invariant mass 
resolution. The signal parameters for these fine mass 
steps are interpolated from the fully simulated samples. 

Figure[3]shows the reconstructed diphoton mass spec- 
trum. No excess is visible. This is quantified by the p- 
value of the background-only hypothesis, which gives 
the fraction of background-only experiments that would 
have a profile likelihood ratio of the zero signal hypoth- 
esis relative to the best-fitted signal strength at least as 
low as the one found in the data. Negative signals are 
not allowed in the fit, so /^-values above 0.5 are trun- 
cated. This /5-value is shown in FigureQa) as a function 
of the hypothesised Higgs boson mass. The minimal p- 
value, corresponding to the largest background upward 
fluctuation, is «5% and is found for a hypothesised mass 
of =s 128 GeV. The probability for such an excess to ap- 
pear anywhere in the investigated mass range is around 
40%. 
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Figure 3: Distribution of the reconstructed diphoton mass. All five 
diphoton categories have been combined. The exponential fit to the 
full sample of the background-only hypothesis, as well as the expected 
signal for a Higgs boson mass of 120 GeV with five times the Standard 
Model predicted yield, are also shown for illustration. 



Exclusion limits on the inclusive production cross- 
section of a Standard Model Higgs boson relative to 
the Standard Model cross-section are derived. For this 
purpose, a modified frequentist approach CL S l52ll . cor- 
responding to the ratio of p-values for the signal-plus- 
background and the background-only hypothesis for a 
given assumed signal strength is used. A given signal 
strength is excluded at 95% confidence level if its CL S 
is smaller than 0.05. The confidence levels are com- 
puted using a large number of signal-plus-background 
pseudo-experiments, with different signal yields, and 
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background-only pseudo-experiments. The results, in- 
cluding systematic uncertainties, are shown in Fig- 
ure Hfb). The impact of the experimental systematic 
uncertainties is about 5% on the expected limit for the 
mass hypothesis of 120 GeV. 

The expected median limit in the case of no signal 
varies from 3.3 to 5.8 as a function of the Higgs boson 
mass. The variations of the observed limit, between 2.0 
and 5.8, are consistent with expected statistical fluctua- 
tions around the median limit. 
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Figure 4: (a) p-value (po) of the background-only hypothesis as a 
function of the investigated Higgs boson mass (b) 95% confidence 
level upper limits on a Standard Model Higgs boson production cross- 
section, relative to the Standard Model cross-section, as a function of 
Higgs boson mass hypothesis. The solid line shows the observed limit. 
The dashed line shows the median expected limit for background-only 
pseudo-experiments, with the bands indicating the expected fluctua- 
tions around this limit at the 1 <x and 2 <x levels. 



7. Conclusions 

A search for the Standard Model Higgs boson in the 
H — > yy decay mode has been performed using an inte- 
grated luminosity of 1.08 fb -1 recorded by the ATLAS 
experiment in 2011. A high purity diphoton sample is 
selected. No excess is found in the diphoton invariant 
mass distribution in the mass range of 110-150 GeV. 



The observed limit on the cross-section of a Standard 
Model-like Higgs boson decaying into a pair of photons 
ranges between 2.0 and 5.8 times the Standard Model 
cross-section. These variations are compatible with sta- 
tistical fluctuations around the expected limit for this 
data set. These limits are the most stringent to date in 
this channel. 
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